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SUMMARY 


A  modular  model  of  the  sudden-expansion  combustor  has  been 
developed  and  applied  to  the  computation  of  a  variety  of 
sudden-expansion  flowfields,  both  reacting  and  non-reacting. 
The  modular  concept  involves  the  representation  of  a  complex 
flowfield  through  the  use  of  detailed  models  of  flowfield 
regions,  coupled  together  through  their  boundary  conditions. 
For  the  sudden-expansion  combustor,  basic  elements  of  the 
model  include  a  parabolic  directed-flow  computation  and  a 
well-stirred  reactor  model  representing  the  recirculation 
region(s).  The  parabolic  module  incorporates  a  two-equation 
turbulent  kinetic  energy  turbulence  model,  and  both  the 
directed  flow  and  stirred  reactor  elements  make  use  of  the 
quasiglobal  formulation  for  the  treatment  of  finite-rate 
hydrocarbon  chemical  kinetics.  Coupling  between  the  modules 
1  is  through  a  shear  layer  representation  which  provides  the 

|  parabolic-flow  boundary  conditions  and  the  stirred  reactor 

'  feed  rates.  Compared  to  unified  elliptic  models  of  the 

1  dump  combustor,  the  major  advantage  of  this  approach  is  the 

greatly  increased  detail  possible  in  the  combustor  analysis. 
This  detail  in  turn  allows  the  development  of  chemical 
i  kinetic  and  flowfield  models  applicable  to  the  refinement 

i  of  elliptic  formulations. 


ABSTRACT 


The  development  and  application  of  a  modular  model  for  the 
prediction  of  the  performance  of  sudden-expansion  burners  as 
a  function  of  the  controllable  parameters  relevant  to  combus¬ 
tor  design  is  described.  The  model  is  based  on  a  concept  in 
which  the  recirculation  zone,  treated  as  a  stirred  reactor, 
is  coupled  to  a  parabolic  boundary  layer  formulation  for  the 
flow  outside  the  recirculation  zone.  Hydrocarbon  oxidation 
kinetics  and  turbulent  kinetic  energy  models  are  employed  in 
the  model  development.  In  addition  to  the  parabolic-flow  and 
stirred  reactor  elements,  a  module  representing  the  fuel  in¬ 
jection  process  has  been  developed.  Results  of  the  applica¬ 
tion  of  the  modular  model  to  the  analysis  of  cold-flow  and 
reacting-flow  dump  combustor  experimental  data  are  described. 
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SECTION  1  -  INTRODUCTION 


Control  of  flame  stabilization  and  flame  propagation  in  a 
turbulent  flow  represents  a  key  element  in  combustion  chamber  design. 
The  placement  and  geometry  of  fuel  injectors,  flameholders ,  and  air 
distribution  ports  are  basic  design  parameters  that  govern  the  per¬ 
formance  of  a  particular  combustor.  Therefore  it  is  desirable  to 
compute  combustion  chamber  flowfields  in  order  to  understand  the 
phenomena  that  occur  in  existing  combustors  and  to  predict  the  per¬ 
formance  of  new  combustor  design  concepts.  In  addition  to  the  over¬ 
all  combustor  flowfield,  modeling  can  also  be  used  to  provide  insight 
into  the  behavior  of  portions  of  the  combustor,  such  as  flameholders 
and  fuel  injectors,  under  a  variety  of  conditions.  The  insight 
gained  through  the  use  of  these  "unit"  analyses  can  be  of  substantial 
use  in  the  planning  of  a  combustor  test  program  and  in  the  interpre¬ 
tation  of  combustor  and  combustor  component  test  data. 

The  computation  of  a  generalized  combustor  flowfield  is  a 
formidable  task,  involving  a  number  of  complex,  coupled  physical  and 
chemical  processes,  which  can  include  turbulent,  recirculating  flow, 
possibly  with  swirl,  finite-rate  chemical  kinetics,  and  droplet 
evaporation  and  combustion.  The  obvious  difficulties  include  not 
only  sorting  out  the  multitude  of  coupled  mechanisms  but  also  involve 
the  typical  disparity  in  characteristic  length  and  time  scales  in 
combustion  chamber  flows.  Despite  the  problems  involved,  considerable 
progress  has  been  made  in  recent  years  in  the  development  of  calcula¬ 
tion  methods  for  these  flows,  using  techniques  for  the  solution  of 
the  elliptic  governing  equations  and  simpler  yet  physically  perceptive 
"modular"  modeling  techniques. 

A  "unified"  model  of  an  overall  combustion  chamber  flowfield 
in  general  requires  the  numerical  solution  of  the  elliptic  form  of  the 
governing  equations,  since  most  practical  combustion  chamber  flows 
involve  large  regions  of  recirculation,  in  which  axial  diffusion  is 
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important.  A  considerable  amount  of  research  effort  has  been  put  into 
the  development  of  numerical  techniques  for  these  problems,  and  success¬ 
ful  comparison  of  calculation  with  experiment  for  recirculating  flows 
with  large  heat  release  has  been  reported  by  Hutchinson,  et  al.  [l]  and 
by  Abou  Ellail,  et  al .  [2],  for  example.  Significantly,  in  both  cases 
a  careful  adaptation  of  the  numerical  model  to  the  specific  experimental 
configuration  was  reported  to  be  required,  and  in  both  of  these  papers 
it  was  noted  that  the  details  of  the  computation  required  careful  han¬ 
dling  to  obtain  the  accuracy  demonstrated.  Furthermore,  Abou  Ellail, 
et  al.  [?]  note  that  it  is  not  possible  to  provide  sufficient  resolution 
in  a  detailed  combustor  flowfield  computation  to  adequately  describe 
processes  such  as  fuel  injection  for  which  the  mixing  process  initially 
occurs  on  a  scale  much  smaller  than  that  of  the  overall  combustion 
chamber. 

Thus,  while  the  development  of  numerical  models  capable  of 
providing  direct  solution  of  the  equations  governing  specific  combustion 
chamber  flowfields  continues,  a  need  exists  for  the  development  of  phys¬ 
ically  perceptive  yet  mathematically  simpler  models.  This  requirement 
arises  from  both  the  need  to  provide  a  model  which  allows  reasonably 
rapid  computation  of  a  number  of  different,  complex  combustor  geometries, 
and  the  need  to  develop  models  for  those  processes,  such  as  fuel  injec¬ 
tion,  which  occur  on  scales  smaller  than  can  be  adequately  resolved  in 
a  detailed  overall  flowfield  computation.  The  development  of  approxi¬ 
mate  methods  -  modular  models  -  is  a  response  to  the  requirements  just 
outl ined. 

The  basic  interest  in  the  application  of  approximate  techniques 
is  to  avoid  the  complexities  inherent  in  a  direct  calculation  of  an 
elliptic  flowfield  by  making  suitable  assumptions  that  allow  the  flow  to 
be  computed  using  simpler  approaches.  Clearly  the  simplest  possible 
procedure  is  to  assume  that  the  flowfield  is  effectively  one-dimensional 
thus  avoiding  any  necessity  for  definition  or  calculation  of  velocity  or 
species  profile  effects.  A  somewhat  more  sophisticated  approach  is  to 
assume  that  the  combustor  flowfield  can  be  broken  down  into  separate 
zones,  each  of  which  can  be  calculated  individually  in  some  detail,  and 


then  coupled  together  in  some  fashion  to  obtain  an  overall  computational 
analog  of  the  combustor  flow.  Such  approaches  are  termed  modular  models, 
examples  of  which  have  been  reported  by  Roberts,  et  al.  [3]  and  Swithen- 
bank,  et  al.  [4];  the  formulation  described  in  this  report,  while  similar 
to  those  of  Roberts,  et  al .  and  Swithenbank,  et  al .  in  that  the  overall 
combustor  flowfield  is  broken  down  into  computational  subunits,  provides 
far  more  detail  than  either  of  the  former  models. 

The  sub-models  used  for  the  flowfield  regions  and  the  assump¬ 
tions  involved  in  the  model  formulations  in  each  of  the  approaches  de¬ 
scribed  in  references  3  and  4  have  been  discussed  by  Harsha  and  Edelman 
[5]-  In  the  present  approach,  the  combustor  flowfield,  represented 
schematically  in  Figure  1,  is  broken  down  into  three  major  components: 
a  directed  flow,  which  is  treated  as  parabolic,  a  recirculation  zone, 
assumed  to  be  represented  by  well-stirred  reactor(s),  and  a  turbulent 
shear  layer  along  the  dividing  streamline  which  separates  the  other  two 
regions.  The  shear  layer  serves  as  the  coupling  region  between  the  other 
two  model  components;  fluxes  of  species  and  energy  across  this  shear 
layer  form  the  boundary  conditions  on  the  two  computational  regions. 
Finite-rate  chemistry,  based  on  the  quasiglobal  model  [6]  is  included  in 
the  formulations  for  both  the  directed  flow  and  well -stirred  reactor  re¬ 
gions,  although  for  the  modular  model  calculations  described  in  this 
report,  the  recirculation  region  well -stirred  reactor  formulation  has 
been  restricted  to  a  global  finite-rate  chemistry  model.  The  directed 
flow  is  assumed  to  be  fully  turbulent,  with  the  turbulent  viscosity  de¬ 
fined  by  a  two-equation  turbulence  model  [7]. 

A  key  difference  between  this  modular  model  and  the  models  of 
Roberts,  et  al .  [3]  and  Swithenbank,  et  al .  [4]  is  the  provision  for  the 
shear  layer  coupling  region  in  the  current  model.  Through  the  use  of 
this  element  of  the  model,  the  division  of  the  mass  flux  between  the 
directed  flow  and  the  recirculation  region  is  computed  iteratively 
rather  than  spen'ficd  empirically.  Furthermore,  the  directed  flow  re¬ 
gion  is  computed  in  detail  as  a  two-dimensional  parabolic  flowfield, 
rather  than  through  a  one-dimensional  approximation,  allowing  the  use 
of  detailed  computations  of  the  mixing  and  chemical  reactions  in  this 
region  of  the  combustor.  • 


Figure  1.  SCHEMATIC  OF  SUDDEN- EXPANSION  (DUMP)  BURNER 

The  features  of  the  three  modular  models  discussed  in  this 
section  are  summarized  in  Table  1.  Although  the  present  model  is  compu¬ 
tationally  more  complex  than  the  models  of  Roberts,  et  al.  [3]  and 
Swithenbank,  et  al.  [4],  and  requires  an  iterative  solution,  the  addi¬ 
tional  complexity  allows  the  development  of  an  overall  combustor  flow- 
field  computation  which  can  include  far  more  detailed  modeling  than  is 
possible  with  either  of  these  approaches.  Furthermore,  each  element  of 
the  present  modular  model  -  the  well -stirred  reactor,  the  parabolic  flow 
computation,  and  the  shear  layer  representation  can  be  systematically 


developed.  For  example,  the  effects  of  unmixedness  can  be  introduced 
into  both  the  stirred  reactor  module  and  the  parabolic  directed  flow 
computation.  Multiple  well-stirred  reactors  can  also  be  easily  intro¬ 
duced  for  more  detailed  modeling  of  the  recirculation  region.  The  detail 
available  in  the  parabolic  flowfield  computation  allows  the  inclusion  of 
models  for  processes  such  as  liquid  fuel  injection,  and  the  spread  and 
vaporization  of  the  resulting  spray,  treated  either  in  bulk  or  through 
the  use  of  detailed  spray  modeling.  Thus,  modular  modeling  of  the  type 
described  in  this  report  offers  considerable  potential  for  systematic 
development  in  complexity  with  the  prospect  of  providing  computational 
capabilities  substantial ly  equivalent  to  those  available  from  direct 
solution  of  the  elliptic  equations,  while  at  the  same  time  providing 
greater  insight  into  the  physical  details  involved  in  a  generalized 
combustor  flowfield. 


Table  1.  COMPARISON  OF  MODULAR  MODELS 


.  ,  r- 

Roterts  ,  el  al .  [ 3j 

Sul teenbanlt,  et  *1 .  [4] 

Present  Model 

Clements 

Plug  Flo*  Reactors 

Plug  F 1  ow.'Pet-  fee  tty 

Stirred  Reactors 

Parabolic  Flow/Sttrred 
Reactors 

Coupling 

Condition, 

empirical  Flow  Division 

Empirical  Flow  Division 

Iterative  Solution  For 
Flow  Division 

Turbulent  flow 
Character! /alien 

None 

None  *  Enters  Empirical 
Flow  Division 

Turbulent  Kinetic 

Energy  Model 

Chemistry  Model 

Quasi  jlohal  M/C 
nineties 

Global  H/C  Kinetics 

Quaslglobal  H/C 

Kinetics 

In  this  report,  the  development  of  the  modular  model  for  a 
sudden-expansion  combustor  is  described  in  detail.  The  model  formula¬ 
tion,  including  the  definition  of  coupling  conditions  and  a  description 
of  the  overall  flowfield  iteration  procedure,  is  described  in  Section  2. 
Results  of  the  application  of  the  approach  to  the  computation  of  a 
variety  of  sudden-expansion  flowfields  are  described  in  Section  3;  where 
the  capabilities  of  this  approach  and  current  unified  elliptic  formula¬ 
tions  overlap,  such  as  in  the  computation  of  incompressible  sudden- 
expansion  flows,  the  results  of  the  modular  approach  are  compared  with 
both  available  experimental  data  and  unified  model  computational  results. 
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As  noted  earlier,  the  analysis  of  combustor  processes  can  make 
use  of  unit  models  as  well  as  the  overall  modular  technique.  These  unit 
models  can  then  in  turn  provide  elements  of  the  modular  approach,  as 
well  as  providing  insight  into  combustor  design  problems  on  their  own. 
Examples  of  the  application  of  unit  process  modeling  are  discussed  in 
Section  4:  these  examples  include  the  use  of  a  simple  fuel  droplet 
vaporization  model  to  interpret  liquid  spray  combustion  efficiency  re¬ 
sults  and  the  use  of  the  well-stirred  reactor  formulation  to  analyze 
flameholder  blowout  problems  and  to  extend  the  range  of  available  blow¬ 
out  correlations  into  regions  applicable  to  ramjet  development  programs. 

Overall  conclusions  reached  from  this  work  with  regard  to  the 
problems  of  combustor  modeling  are  discussed  in  Section  5.  The  work 
described  in  this  report  shows  that  there  are  considerable  advantages  to 
be  gained  from  the  use  of  analytical  models  as  interpretative  techniques 
and  planning  guides  in  conjunction  with  experimental  combustor  develop¬ 
ment  programs,  such  as  the  coordination  between  this  work  and  the  experi¬ 
mental  work  at  the  Air  Force  Aeropropulsion  Laboratories  which  has  been 
maintained  throughout  the  effort  described  herein. 


SECTION  2  -  MODULAR  MODEL  FORMULATION 

The  basic  elements  of  the  current  modular  model  formulation 
are  the  parabolic  finite-difference  computational  technique  developed 
by  Boccio,  Weilerstein  and  Edelman  [8],  modified  to  incorporate  the  two- 
equation  turbulent  kinetic  energy  model  developed  by  Launder,  et  al .  [7], 
used  for  the  directed-flow  portion  of  the  analysis,  and  the  stirred  re¬ 
actor  computation  developed  by  Edelman  and  Weilerstein  [9].  Both  of 
these  elements  make  use  of  the  quasiglobal  model  developed  by  Edelman 
and  Fortune  [6]  for  rapid  computation  of  finite-rate  hydrocarbon  air 
kinetics.  These  elements,  or  modules,  are  coupled  together  through  a 
simplified  representation  of  the  turbulent  shear  layer  which  exists  be¬ 
tween  the  directed  flow  and  the  recirculation  region.  In  the  modular 
approach,  the  shear  layer  representation  is  used  to  define  the  gradients 
in  velocity,  species,  and  enthalpy  between  the  two  regions  of  the  flow, 
thus  providing  both  the  boundary  conditions  on  the  directed  flow  and  the 
stirred  reactor  feed  rates.  Oetails  of  the  mathematical  formulations 
used  for  the  flowfield  regions  are  discussed  in  this  section;  details  of 
the  numerical  analysis  procedure  may  be  found  in  Appendix  I. 

2.1  WELL-STIRRED  REACTOR  -  THE  RECIRCULATION  ZONE  MODEL 

Flowfield  regions  in  which  intense  backmixing  occurs  can  ap¬ 
proach  the  limit  of  complete  mixing,  and  thus  the  well -stirred  reactor 
concept  is  attractive  for  representing  the  recirculation  region(s)  in  a 
sudden-expansion  burner.  For  'a  well-stirred  reactor,  the  equations  de¬ 
scribing  the  transport  of  energy  and  species  reduce  to  the  relations 

Continuity 

m*  =  m®  =  m  .  (1) 
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where  the  superscripts  I  and  0  represent  inflow  to  and  outflow  from  the 

stirred  reactor,  M  the  mass  contained  within  the  stirred  reactor,  a.  are 

.  t  1 

the  mass  fractions  of  the  i  species  in  the  reactor,  the  net  inflow  of 

species  i,  V  is  the  reactor  volume,  w.  the  volumetric  production  rate  for 

*  • 

species  i  due  to  chemical  reaction,  h  the  enthalpy,  and  Q  the  rate  of  heat 
addition  i.o  the  reactor  from  its  surroundings. 

Note  that  in  this  set  of  equations  the  species  transport  equa¬ 
tion  is  written  in  non-steady  form.  This  formulation  has  been  adopted 
to  facilitate  solution  of  the  stirred  reactor  governing  equations  with 
finite-rate  chemical  kinetics;  the  steady-state  stirred  reactor  solution 
is  obtained  when 


d  a. 
_ 

dt 


=  0 


(4) 


These  equations,  along  with  expressions  for  the  volumetric  pro¬ 
duction  rates  appearing  in  Eq:  2,  expressions  for  the  enthalpy  as  a 
function  of  species  concentrations,  and  the  equation  of  state,  define  the 
temperature  and  species  concentrations  in  the  stirred  reactor,  given  the 
inflow  rates  for  species  and  enthalpy.  In  the  modular  model,  the  net 
inflow  of  species  and  enthalpy  can  be  expressed  as  line  integrals  involving 
gradients  evaluated  along  the  dividing  streamline,  so  that  for  the  modular 
model  the  energy  and  species  conservation  equations  for  the  stirred  reactor 
can  be  written 


where  Rc(x)  represents  the  radial  location  of  the  dividing  streamline, 
and 


Here  the  superscript  I  refers  to  inflow  into  the  stirred  reactor  (re¬ 
circulation  region)  and  superscript  0  refers  to  outflow  from  the  stirred 
reactor,  V  is  the  reactor  volume,  pc  a  characteristic  density  of  the 
stirred  reactor  region,  and  pv-j.  a  characteristic  eddy  diffusivity,  evalu¬ 
ated  from  the  outer  flowfield  solution  in  the  region  of  the  dividing 

streamline.  The  term  w-  represents  the  rate  of  production  of  species  i 

1  « 

caused  by  chemical  reactions  and  Q  represents  the  heat  input  to  the 
stirred  reactor  region  through  the  combustor  walls. 

2.2  PARABOLIC  MIXING:  THE  DIRECTED-FLOW  MODEL 

The  second  major  element  of  the  modular  model  for  a  sudden- 
expansion  combustor  is  the  formulation  for  the  directed  flow  portion  of 
the  combustor  flowfield.  It  is  assumed  that  the  boundary  layer  approxi¬ 
mations  apply  to  this  part  of  the  flowfield,  so  that  the  des.cribing 
equations  are  parabolic.  For  a  steady,  axisymmetric  flow  these  equations 
may  be  written: 
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Global  Co n tinui ty 


3rpu  cirp^  _  n 

3x  ~3r  u 


Species  Diffusion  for  the  ith  Specie 


3ot.  3a-  i  3 
pu  TT  +  pv  TF  =  r  3r  |rp  Sc  ' 


Mb, 


Momentum  Equation : 


3u  ,  3u  _  1  (  3  /  3u  \  1  3p 
PU3^+  pV3?~  Fl  3F(rpvT  r  "  9X 


(7) 


(8) 


(9) 


and  the  Energy  Equation 


pu 


3H 

3x 


in  which  Pr  and  Sc  represent  the  Prandtl  and  Schmidt  numbers,  respectively. 
These  equations,  along  with  expressions  for  the  enthalpy 


H  =  h  +  Hj-  and  h  =  £  h^  (T) 


(11) 


and  the  equation  of  state 


p  *  pRT  £(a./W.)  (12) 

i 

can  be  solved,  given  an  expression  for  the  turbulent  momentum  diffusivity 
Vj,  or  the  turbulent  eddy  viscosity  p-j-  =  pvy.  The  problem  of  the  most 
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appropriate  formulation  for  the  eddy  viscosity  in  a  turbulent  flow  has 
occupied  the  attention  of  researchers  in  turbulent  flow  for  many  years, 
and  numerous  proposals  for  the  appropriate  form  of  the  eddy  viscosity 
have  been  made.  In  general,  up  until  a  few  years  ago,  models  for  the 
turbulent  eddy  viscosity  involved  a  relation  between  a  local  length  scale 
and  a  local  measure  of  the  velocity  gradient.  The  free  turbulent  mixing 
model  proposed  by  Prandtl  in  1942  [lO]  is  a  case  in  point;  with  this 
model 


Pj  =  pcb | Au |  (13) 

where  c  is  a  constant,  b  is  a  measure  of  the  width  of  the  mixing  region, 
and  Au  a  measure  of  the  velocity  difference  across  the  mixing  region. 

While  eddy  viscosity  models  have  in  certain  circumstances 
enabled  successful  calculations  of  particular  turbulent  flows  to  be  car¬ 
ried  out,  in  general,  the  constant  of  proportionality  involved  in  the 
model  varies  markedly  and  unpredictably  in  different  flowfields  [l 1 ]. 

The  observed  lack  of  generality  of  eddy  viscosity  models  has 
led  to  the  development  of  a  class  of  models  in  which  additional  partial 
differential  equations  are  written  to  obtain  the  spatial  variation  of 
the  shear  stress.  Of  these,  one  of  the  most  highly  developed  is  the 
two-equation  turbulent  kinetic  energy  model  developed  by  Launder,  et  al . 
[7],  in  which  transport  equations  are  written  for  the  turbulent  kinetic 
energy 

k  =  j  (u'2  +  v'2  +  w'2)  (14) 

and  its  dissipation  rate,  e.  In  boundary  layer  form,  these  equations  can 
be  written: 

Turbulent  Kinetic  Energy 
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Turbulence  Energy  Dj ssipation 


where 


pu 


3  e  , 

3  c 

1 

3  1 

(bT»' 

"37  + 

pv  w  = 

r 

Dr  1 

l  a 
\  e 

CE1  kMT 


(16) 


bj  =  p k2/e 


(17) 


and  in  which 


=  0.09,  CE1  =  1.40,  CE2  =  1.95  (18) 

The  turbulence  energy  dissipation  rate,  e,  can  be  related  to  a  turbulent 
length  scale,  2.^,  through  the  Kolmogorov  hypothesis 

e  oe  k3/2/£k  (19) 

so  that  the  transport  equation  for  turbulence  energy  dissipation  can  also 
be  regarded  as  a  description  of  the  spatial  variation  of  the  turbulent 
length  scale. 

2.3  CHEMICAL  KINETICS:  THE  QUASIGLOBAL  MODEL 

In  both  the  stirred  reactor  and  directed  flow  portions  of  the 
modular  model  volumetric  chemical  kinetic  production  rate  terms  appear 
in  the  species  transport  equations.  These  terms  can  be  evaluated  using 
a  full  hydrocarbon  chemical  kinetics  scheme  based  on  the  quasiglobal 
kinetics  model  [6],  which  has  as  a  key  element  a  subglobal  oxidation  step 

CnHm  +  7  °2  +  I  H2  +  nC0  (20) 

This  reaction  is  unidirectional  with  an  empirically  determined  rate  (grams 
of  fuel/cc/sec)  given  by 
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ATbP°'3  [CnHni]  "  [02]  exp  [-  E/RT  ]  (21) 

with  the  constants  A,  b,  and  E/R  defined  in  Table  ?,  where  P  must  be  given 
in  atmospheres,  T  in  degrees  Kelvin  and  [  ]  denotes  molar  concentration; 
coupled  to  this  subglobal  step  are  the  intermediate  reversible  reactions 
given  in  Table  2. 

Table  2.  EXTENDED  C-H-0  CHEMICAL  KINETIC  REACTION  MECHANISM 

kf  =  ATbexp  (-E/RT) 


fraction 

A 

forward 

E/P 

D 

wh 

*  Jyn  CO* 

A.  Long  Chain  6.0  x  10*- 

8.  Cyclic  2.08  x  10' 

1 

1 

12.2  x  10, 
19.65  x  10J 

2) 

CO  4  OH 

-  H  ♦  C02 

5.6  x  1011 

0 

0.543  x  103 

3) 

CO  ♦  o2 

-  co2  ♦  0 

3.0  x  )012 

0 

25.0  x  103 

4) 

CO  ♦  0  4 

H  •  C02  ^  M 

1.8  x  1019 

-1 

2.0  x  10J 

5) 

k2  ♦  02 

•  OH  ♦  OH 

1.7  x  1013 

0 

24.7  x  10J 

6) 

OH  4  H2 

-  H20  ♦  H 

2.19  x  10U 

0 

2.59  x  105 

7) 

OH  4  OH 

-  0  ♦  H2O 

5.75  x  1012 

0 

0.393  x  103 

8) 

0  4  H2  . 

H  4  OH 

1.74  X  1013 

0 

4.75  x  103 

9) 

H  i  0;  ■ 

0  4  OH 

.  2.24  x  101* 

0 

8.45  x  103 

10)  H  4  0  4  H  -  OH  ♦  M 

1.0  x  )DU 

0 

D 

11)  H  ♦  0  ♦  0  •  02  ♦  H 

9.38  x  101* 

0 

0 

12)  N  4  H  4  H  ■  H2  4  H 

6.0  x  1015 

0 

0 

13)  H  *  H  4  OH  ■  H20  4  « 

1.0  x  1017 

0 

0 

14)  0  ♦  Nj  -  K  4  NO 

1.36  x  !0U 

0 

3.775  x  10* 

15)  N?  4  02  •  N  4  N02 

2.7  x  101* 

-i.e 

6.06  x  10* 

16)  N2  4  02  •  NO  ♦  NO 

9.1  x  102* 

-2.5 

6.46  x  10* 

17)  NO  •»  NO  •  N  4  NOj 

1.0  x  tO10 

0 

4.43  x  10* 

18)  NO  4  0  -  02  4  N 

1.55  x  109 

1.0 

1.945  x  10* 

19)  N  4  NO  •  0  4  N  4  H 

2.27  x  1017 

-0.5 

7.49  i  10* 

20)  M  ♦  N0?  -  0  ♦  NO  ♦  H 

1.1  x  10U 

0 

3.30  I  10* 

21)  K  4  N02  •  02  *  N  a  W 

6.0  x  101* 

-1.5 

5.28  x  10* 

22)  NO  ♦  02  -  N02  *  0 

1.0  x  1012 

0 

2.29  X  10* 

23)  N  4  OH  •  NO  4  H 

4.0  x  1013 

0 

0 

24)  N  4  NOj  •  NO  +  OH 

3.0  i  1013 

0 

0 

25)  C02  4  N  •  CO  ♦  NO 

.  2.0  x  1011 

-1/2 

4.0  x  103 

26)  CO  4  NOj  •  C02  4  NO 

2.0  x  1011 

-1/2 

2.5  x  103 

[c]  •  [,)  -°k.  tp]  - to  •  Si 

reaction  rate  It  obtained  from  kf  md  the  equilibrium  constant  K* 
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2.4 


MODELS  FOR  THE  FUEL  INJECTION  PROCESS 


In  addition  to  the  basic  components  of  a  parabolic,  directed 
flow  analysis  and  a  well-stirred  reactor  formulation,  the  modular  concept 
can  be  extended  to  include  modules  which  represent  other  elements  of  the 
dump  combustor  flowfield,  for  example,  the  fuel  injection  process.  The 
detail  of  the  computation  provided  by  the  use  of  a  parabolic  directed  flow 
analysis  (as  opposed  to  the  relative  coarseness  of  the  numerical  grid 
allowable  in  current  unified,  elliptic  solution  techniques)  is  the  key 
feature  of  the  modular  model  that  allows  the  inclusion  of  a  fuel  injection 
module  in  the  complete  analysis.  This  is  particularly  true  in  the  case 
of  liquid  fuel  injection,  for  at  the  fuel/air  ratios  appropriate  for  dump 
combustor  operation,  the  liquid  fuel  streams  initially  occupy  a  very 
small  portion  of  the  overall  combustor  cross-sectional  area.  It  is  also 
worth  noting  that  it  has  been  stated  that  the  further  development  of  uni¬ 
fied,  elliptic  combustor  flowfield  calculations  requires  the  development 
of  what  are  essentially  modular  models  for  such  features  as  the  fuel  in¬ 
jection  process  [2]  since  these  key  processes  occur  on  scales  much  smaller 
than  feasible  numerical  resolution  allows.  Thus  the  development  of  fuel 
injection,  vaporization,  and  spreading  models  through  use  of  the  overal 
modular  model  concept  is  necessary  for  further  development  of  unified 
techniques. 

The  liquid  fuel  injection  model  makes  use  of  a  combination  of 
empirical  information  and  turbulent  mixing  calculations.  For  example, 
the  fuel  jet  penetration  from  the  wall  is  computed  through  the  use  of  the 
penetration  correlation  developed  by  Catton,  Hill  and  McRae  [12],  using 
the  breakup  time  correlation  developed  by  Clark  [ 1 3]  to  compute  the  down¬ 
stream  distance  at  which  penetration  is  to  be  computed.  That  is,  it  is 
assumed  that  the  fuel  jet  has  turned  parallel  to  the  airflow  at  the  axial 
position  at  which  the  initial  fuel  jet  has  broken  up  into  droplets,  as 
given  by  the  breakup  time  correlation  and  the  local  airflow  velocity. 

Since  the  basic  modular  model  formulation  involves  an  axisym- 
metric  flowfield,  individual  fuel  jets  cannot  be  resolved,  and  it  is 
assumed  in  the  formulation  that  the  liquid  fuel  spray  forms  an  annulus 
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whose  cross-sectional  area  may  either  be  specified  or  computed  based  on 
an  assumed  fuel  spray  bulk  velocity.  A  bulk  spray  evaporation  correla¬ 
tion  is  then  used  to  compute  the  fuel  vaporization  rate;  this  correlation, 
developed  by  Ingebo  and  Foster  [14]  is  a  function  of  the  initial  velocity 
and  temperature  difference  between  the  fuel  spray  and  the  surrounding  air 
stream.  Spreading  of  the  fuel  jet  is  computed  through  use  of  a  turbulent 
mixing  hypothesis  as  for  the  mixing  process  in  the  remainder  of  the  para¬ 
bolic  flow. 

Figure  2  shows  the  results  of  a  computation  of  the  fuel  injec¬ 
tion  process  for  three  fuel  injectors,  located  in  the  combustor  inlet 
wall,  along  the  centerline,  and  in  a  midstream  position.  Shown  are  the 
computed  contours  of  the  fuel  mass  fraction,  cip,  with  the  vapor-phase 
fuel  shown  as  the  solid  line  and  the  liquid  phase  fuel  as  the  dotted  line, 
as  a  function  of  both  axial  and  radial  position  in  the  combustor  inlet. 

For  these  calculations  the  air  inlet  velocity  was  approximately  700  ft/sec 
at  a  temperature  of  1600°K;  the  overall  fuel-air  equivalence  ratio  was 
0.6.  The  results  shown  in  Figure  2  provide  a  good  example  of  the  detail 
of  the  fuel  injection  process  available  through  use  of  this  aspect  of 
modular  modeling. 

2.5  COUPLING  RELATIONS  AND  ITERATION  PROCEDURE 

The  overall  logic  of  the  calculation  procedure  is  shown  schemat¬ 
ically  in  Figure  3.  Inlet  conditions  can  be  defined  upstream  or  down¬ 
stream  of  tiie  fuel  injection  station  as  required;  if  liquid  fuel  injection 
is  specified,  the  penetration  and  spreading  of  the  liquid  fuel  jets  are 
computed  as  described  in  the  preceding  section.  The  computation  proceeds 
to  the  dump  plane  in  the  sudden  expansion  combustor,  at  which  point  the 
dependent  variables  in  the  flowfield  are  stored  for  later  use  in  the  re¬ 
circulation  zone  iteration  procedure. 

An  initial  state  for  the  stirred  reactor  must  be  specified  in 
order  to  begin  the  iteration  procedure.  This  state  is  reasonably  arbi¬ 
trary,  except  that  a  reacted  temperature  level  must  be  specified.  The 
shape  of  the  dividing  streamline  separating  the  directed  flow  and  the 


-15- 


RADIAL  POSITION,  FT 


.6 


.5  -  0  a  -5  0 - 75 

a 


r  f 

- 1  , 

n  t;  - 1 


0  a  -5  0  .5  L - J 

a  0  .5 

a 


INLET  AIR  TEMPERATURE  =  1600  K 

INLET  AIR  VELOCITY  =  700  ft/sec  _ ^ 

OVERALL  EQUIVALENCE  RATIO  =  0.6 


^LIQUID  PHASE 
-VAPOR  PHASE 


AXIAL  POSITION,  FT. 


Figure  2.  PREDICTED  FUEL  MASS  FRACTION  DISTRIBUTIONS  FOR 
SIMULTANEOUS  INJECTION  AT  THREE  RADIAL  LOCATIONS 
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recirculation  zone  must  also  be  specified,  along  with  the  shear  stress 
level,  expressed  as  a  "skin  friction"  coefficient.  The  initial  stirred 
reactor  state,  in  conjunction  with  a  model  for  the  shear  layer  existing 
in  the  combustor  flowfield  between  the  recirculation  zone  and  the  directed 
flow  then  defines  the  boundary  conditions  for  the  parabolic  directed  flow 
calculation.  Directed  flow  computations  are  carried  out  to  the  axial  sta¬ 
tion  at  which  the  end  of  the  recirculation  zone  has  been  defined;  as  part 
of  these  computations  the  diffusive  flux  of  species  and  energy  across  the 
dividing  streamline  is  computed.  These  fluxes  then  define  a  new  set  of 
stirred  reactor  "feed  rates",  i.e.,  species  and  energy  fluxes  into  the 
recirculation  zone  region. 


Figure  3.  MODULAR  MODEL  CALCULATION  LOGIC 
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Convergence  of  the  procedure  is  defined  by  the  change  in  the 
stirred  reactor  feed  rates  from  iteration  to  iteration;  each  species  and 
the  energy  flux  must  change  less  than  0.1%  before  the  convergence  criteria 
are  satisfied.  If  they  are  not,  the  species  and  energy  fluxes,  and  the 
overall  diffusive  mass  flux  and  physical  recirculation  zone  volume  are 
used  to  compute  a  new  stirred  reactor  state.  The  pressure  required  for 
this  computation  is  taken  to  be  the  arithmetic  average  of  the  pressures 
computed  as  part  of  the  directed  flow  solution  at  the  beginning  and  end 
of  recirculation  region.  The  new  stirred  reactor  state  is  then  used  to 
define  new  directed  flow  boundary  conditions,  and  the  computation  is  re¬ 
started  from  the  dump  station.  When  convergence  is  achieved,  typically 
in  3-5  iterations,  the  combustor  calculation  is  continued  to  the  specified 
combustor  exit  station. 

As  will  be  noted  in  the  description  which  follows,  there  are  a 
number  of  specifiable  constants  in  this  flowfield  model,  including  the 
length  of  the  recirculation  zone,  the  recirculation  zone  shape,  the  initial 
thickness  and  spread  rate  of  the  shear  layer  along  the  dividing  streamline, 
and  the  dividing  streamline  shear  stress  distribution.  Both  the  size  and 
shape  of  the  recirculation  zone  can  be  adequately  specified  as  functions 
of  the  dump  combustor  geometry,  and  it  has  been  found  that  the  shear  layer 
initial  thickness,  spread  rate,  and  shear  stress  distribution  along  the 
dividing  streamline  have  a  marked  effect  on  the  static  pressure  distribu¬ 
tion  computed  for  the  combustor.  This  suggests  that  an  outer  iteration 
loop  can  be  used  to  adjust  the  computation  for  a  specific  combustor  geom¬ 
etry  through  matching  of  the  computed  static  pressure  profile  with  experi¬ 
mental  data  for  that  geometry,  and  as  will  be  described  in  Section  3,  this 
technique  has  been  successfully  used  to  match  computed  and  experimental 
distributions  of  static  pressure  for  both  cold  flow  and  reacting  flow, 
through  adjustment  of  the  specifiable  constants  in  the  model. 

In  describing  the  coupling  between  the  stirred  reactor  and 
directed  flow  elements  of  the  module,  it  is  useful,  for  clarity,  to  con¬ 
sider  a  simplified  formulation  for  the  stirred  reactor,  using  a  global 
finite  rate  kinetics  model.  In  this  model,  the  fuel  consumption  rate  is 
given  by  the  expression 


(22) 


dCr 

■df  =  -CF  C023/2  A  exp  (-B/RT) 

where  the  C's  have  the  units  (moles/cc/sec) ,  T  is  in  °K,  and  R  is  in 
kcal/molc  °K;  A  =  3.80  x  1015  and  B  =  15,600.  Performing'  the  manipula¬ 
tions  necessary  to  recast  Eq.  22  in  terms  of  mass  fractions,  and  substi¬ 
tuting  into  Eq.  2,  results  in  a  set  of  simultaneous  equations  for  the 
active  species  of  the  form: 


dM  a, 


■dt^  =  "‘o2  -  vfV/2''£,<t) 


(23) 


dM  otr  T  aFao//2 

Tt-  -  %  -  WF  TnT^lT  V  6  (T) 


(24) 


dM  acQ 

dt  =  mco2  “  wco2  (r+nsi74')  aFao23/2  V  0  ^ 


(25) 


dM  a, 


Ho0 


dt  =  mH20  +  WH20  ^n  aF  a023/2  V  8  ^ 


(26) 


dM  a. 


2  _  -I 


dt  =  mN 


(27) 


where  n  and  m  are  the  carbon  and  hydrogen  atoms  in  the  general  fuel  mole¬ 
cule  CnHm,  the  W's  are  the  molecular  weights  of  the  species,  and 

6  (T)  =  p3/z  A  exp  (-B/RT)/W  FUn  3/2  (28) 

2 

Given  the  specification  of  the  m! ,  Eq.  23-27,  along  with  the  perfect  gas 
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equation  of  state  and  Eq.  3  are  solved  to  obtain  the  new  steady-state 
values  of  the  active  species,  i.e.,  the  values  obtained  when 

dMcx. 

"I  n 


A  key  element  in  the  modular  formulation  is  the  model  for  the 
shear  layer  between  the  recirculation  zone  and  the  directed  flow.  It  is 
the  computation  in  this  region  that  provides  the  boundary  conditions  for 
the  parabolic  directed  flow  calculation;  through  the  gradients  in  energy 
and  species  mass  fractions  computed  from  the  directed  flow  solution,  the 
influx  rates  in  j  in  Eq.  23-27  and  the  energy  flux  h*  in  Eq.  3  are  evaluated. 

In  the  modular  model  the  shear  layer  region  is  computed  in  a 
simplified  form:  it  is  assumed  that  the  shear  layer  can  be  modeled  as  a 
region  £(x)  across  which  all  dependent  variables  (i.e.,  velocity,  temper¬ 
ature,  species  mass  fractions)  vary  linearly.  It  is  further  assumed  that 
the  shear  layer  width  itself  is  specifiable  by  a  linear  growth  law, 

i  ~  ax  +  b  (29) 

in  which  a  and  b  are  constants.  The  constant  b  can  be  related  to  the 
initial  boundary  layer  thickness  at  the  expansion  plane,  while  a  can  be 
related  to  empirical  expressions  for  shear  layer  growth  rate.  It  might 
further  be  noted  that  for  self-similar  shear  layers,  the  constant  a  can 
be  related  to  the  shear  stress  level  at  some  characteristic  point  in  the 
shear  layer.  However,  for  reacting  sudden-expansion  burner  flows  self¬ 
similar  shear  layer  concepts  do  not  apply  and  the  dividing  streamline 
shear  stress  is,  in  the  modular  model,  independently  specifiable. 

As  shown  in  Figure  4,  the  dividing  streamline  shape  Rc(x)  is 
specified:  this  specification  also  defines  the  dividing  streamline  stream 
function  value  4>w.  Along  \p  =  the  turbulent  shear  stress  is  specified 
through  the  use  of  an  input  "skin  friction  coefficient": 

Cpyj  =  T^/^pU2  (30) 
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where  p  and  D  are  evaluated  as  the  average  of  their  values  along  the  flow- 
field  centerline  and  along  the  streamline  ^  -  iJj  .  As  noted  above,  temper¬ 
ature  and  species  mass  fractions  are  assumed  to  vary  linearly  across  the 
shear  layer.  If  ^  is  defined  as  the  streamline  immediately  "outside"  the 
recirculation  region  (i.e.,  within  the  directed  flow)  and  is  a  (an 
imaginary)  streamline  immediately  "inside"  the  recirculation  zone,  then, 
for  species  i 


aip  ~  air 
3r  /w  “  «Tx) 


(31) 


and  for  temperature 


(32) 


where  a^r  and  Tr  are  the  values  of  the  species  mass  fractions  and  temper¬ 
ature  obtained  from  the  stirred  reactor  solution. 


Figure  4.  DIVIDING  STREAMLINE  DEFINITION  SKETCH 


The  boundary  conditions  at  i|»  =  \plt  established  by  Eq.  30-32 

w 

define,  through  the  parabolic  solution,  effective  gradients  which,  through 
turbulent  diffusion,  can  be  used  to  define  the  influx  into  the  stirred  re¬ 
actor.  For  example,  we  can  write 
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Da, 

ar  Rc(x)dx 


where  L  is  the  distance  measured  along  the  dividing  streamline  and  Vy 
the  turbulent  diffusivity  for  momentum;  Sc  is  the  Schmidt  number.  In 
discretized  form  (Appendix  I),  using  the  transformation 


purdr  =  ip  dtp 


(34) 


Eq.  33  can  be  rewritten 


m: 


I  _  -tt 
Sc 


a?)  Rc(x)  dx 


where  d  n  pur2py/i|j  and  yy  =  pVy.  Defining 


and  noting  that,  by  definition,  aP  is  constant. 
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(37) 


Substituting  Eq.  37  into  Eq.  2,  and  taking  dMo^/dt  =  0  results  in  the 
working  form  of  the  well-stirred  reactor  species  equation  used  in  the 
modular  model , 


•J 


n  +  1 


+  wi  V/ md 


(38) 


where  n  and  n  +  1  are  sequential  elements  of  the  iteration  loop  shown  in 
Figure  3. 
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A  similar  development  can  be  followed  for  the  energy  flux  term 
in  Eq.  3.  For  unity  Lewis  number,  the  result  is 


where  the  summation  is  taken  only  over  influx  species,  i.e.,  species  for 

r  D 
which  a.j  <  ctj  . 

With  the  describing  equations  derived,  the  iteration  procedure 
involved  in  the  modular  model  can  be  restated.  Given  the  conditions  in 
the  directed  flow  at  the  start  of  the  recirculation  zone,  either  by  cal¬ 
culation  of  the  process  occuring  upstream  of  the  dump  station  or  as  input 
conditions,  and  given  initial  guesses  for  the  recirculation  zone  state, 
aC  and  Tr,  as  well  as  the  geometric  variables  Rc(x)  and  a  and  b  in  the 
shear  layer  width  expression,  Eq.  29,  and  the  "skin  friction  coefficient" 
CpW,  the  parabolic  calculation  proceeds  to  the  end  of  the  recirculation 
zone,  using  Eq.  30-32  as  boundary  conditions.  When  the  end  of  the  recir¬ 
culation  zone  is  reached,  aP,  hp,  and  are  evaluated,  and  these  quan¬ 
tities  are  used,  along  with  an  input  specification  of  Q,  to  compute  a 
new  stirred  reactor  state  from  Eq.  23-28.  The  process  is  repeated  until 
changes  smaller  than  0.10%  are  observed  in  aP  and  hp,  after  which  the 
remainder  of  the  combustor  calculation  is  carried  out. 

Boundary  conditions  for  the  turbulent  kinetic  energy,  k,  and 
the  turbulent  kinetic  energy  dissipation  rate,  e,  are  established  in  the 
same  manner  as  described  for  the  mean  total  energy  and  species  equations. 
Thus,  values  are  assumed  for  the  turbulent  kinetic  energy  and  dissipation 
length  scale  for  the  stirred  reactor  region,  kp  and  eR,  respectively,  and 
the  gradient  in  k  and  r  along  the  dividing  streamline  established  from 
the  relations 


9k  \  kP  '  kr 

ar  ;w  ’ 


(40) 


and 
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(41) 


3e\  _  CF  "  er 
3r/w  "  £JxT~ 

As  is  evident  from  this  discussion,  a  number  of  parameters  must 
be  specified,  either  empirically  or  through  comparison  of  model  predic¬ 
tions  with  experiment,  for  a  modular  model  computation  to  proceed.  These 
include  the  size  and  shape  of  the  recirculation  zone,  the  coefficients  a 
and  b  in  Eq.  29,  the  recirculation  zone  values  of  k  and  e,  and  the  divid¬ 
ing  streamline  skin  friction  coefficient  Cp^.  While  this  is  a  fairly 
lengthy  list  of  coefficients,  and  is  in  addition  to  the  inlet  and  wall 
boundary  conditions  that  must  be  specified  in  any  formulation,  experience 
has  shown  that  good  results  can  be  achieved  using  a  reasonably  limited 
range  of  values  of  these  coefficients.  However,  it  should  be  stressed 
that  the  modular  approach  is  not  a  predictive  tool  in  the  sense  that  it 
can  be  relied  upon  in  the  absence  of  experimental  data;  instead,  the 
model  is  one  which,  given  a  limited  amount  of  available  data,  can  be  used 
to  interpret  the  phenomena  occuring  within  the  combustor  under  test,  and 
to  provide  a  means  for  scaling  limited  test  results  to  cover  a  more 
general  range  of  interest  to  the  combustor  designer. 
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SECTION  3 

APPLICATION  OF  MODULAR  MODELING  TO  SUDDEN-EXPANSION  FLOWFIELDS 


In  order  to  demonstrate  the  application  of  the  modular  model  to 
the  analysis  and  interpretation  of  sudden-expansion  flowfield  data,  calcu¬ 
lations  of  a  variety  of  different  sudden-expansion  flows  have  been  carried 
out.  These  flowfields  include  incompressible  and  compressible  nonreacting 
flows  and  premixed  reacting  flows;  in  the  latter  case,  both  one-step 
global  and  the  full  quasiglobal  chemical  kinetics  models  have  been  used. 
Although  there  are  a  number  of  empirical  coefficients  in  the  model,  these 
computations  illustrate  that  coefficient  values  can  be  selected  based  on 
one  set  of  data  and,  in  large  part,  used  for  the  computation  of  different 
flowfields:  where  changes  in  coefficients  are  required  the  differences 
in  each  case  can  be  related  to  physical  differences  in  the  flowfields  in 
different  configurations. 

3.1  COLD  FLOW  APPLICATIONS 

Validation  and  demonstration  of  the  modular  model  requires  de¬ 
tailed  data  for  the  sudden-expansion  configuration  of  interest.  One  such 
source  of  detailed  data  is  the  work  reported  by  Chaturvedi  [15]  for  an 
essentially  incompressible  sudden-expansion  flowfield. 

The  sudden-expansion  configuration  investigated  by  Chaturvedi 
involved  an  area  ratio  of  4.0  with  an  inlet  Reynolds  number  (based  on  in¬ 
let  diameter)  of  2  x  10".  The  inlet  velocity  was  thus  of  the  order  of 
100  ft/sec.  A  variety  of  expansion  angles  were  studied,  varying  from  15° 
to  90°;  the  90°  expansion  angle  case  is  of  the  most  interest  in  this  model 
validation  work.  Chaturvedi 's  apparatus  had  an  L/D  of  25  based  on  the  in¬ 
let  diameter  and  wa s  aspirated  by  a  centrifugal  pump  and  butterfly  valve 
assembly  at  the  downstream  end;  the  inlet  was  a  bellmouth  intake  directly 
coupled  to  the  expansion  section.  Because  of  the  differences  in  Chatur- 
vedi's  experimental  apparatus  relative  to  a  typical  sudden-expansion 
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combustor  configuration,  which  include  the  lack  of  an  exhaust  nozzle 
section,  the  greater  mixing  chamber  length  relative  to  a  typical  combus¬ 
tion  chamber,  and  the  lack  of  an  inlet  section,  the  modeling  assumptions 
necessary  to  match  the  modular  model  predictions  with  these  data  are  not 
necessarily  appropriate  for  a  sudden-expansion  combustor.  Nevertheless, 
Chaturvedi's  data,  which  include  axial  and  radial  profiles  of  mean  veloc¬ 
ity,  axial  and  radial  turbulent  intensity  components,  and  total  pressure, 
are  detailed  enough  to  provide  a  critical  test  of  modeling  capability. 

The  data  used  for  the  comparisons  shown  here  have  been  obtained 
from  that  presented  in  Chaturvedi's  paper  [ 1 5 ] .  Further,  both  the  turbu¬ 
lent  kinetic  energy  and  wall  static  pressure  data  have  been  obtained  by 
calculation  from  the  data  presented  by  Chaturvedi ,  and  are  thus  based  on 
assumptions  not  inherent  in  Chaturvedi's  work.  The  turbulent  kinetic 
energy  has  been  obtained  from  Chaturvedi's  u’/Uq  and  v'/Uq  data  using  the 
assumption  w’/Uq  =  \  v'/Uq,  while  the  wall  static  pressure  has  been  de¬ 
rived  from  measured  centerline  total  pressure  and  velocity  data  using 
Bernoulli's  equation  and  the  assumption  that  the  static  pressure  is  radi¬ 
ally  constant. 

The  parameters,  in  addition  to  wall  boundary  conditions  and  in¬ 
let  conditions,  that  must  be  specified  for  a  modular  computation  include: 

•  recirculation  zone  size  and  shape 

•  shear  layer  initial  thickness  and  growth  rate 

•  dividing  streamline  shear  stress  coefficient 

0  recirculation  zone  turbulent  kinetic  energy  level 

t  recirculation  zone  turbulent  dissipation  length  scale 

Previous  work  had  established  that  the  dividing  streamline  shape  was  well 
represented  by  a  parabolic  arc  arranged  such  that  at  the  step,  dr/dx  =  0, 
and  which  intercepts  the  wall  at  a  distance  9.2  step  heights  from  the  in¬ 
let.  The  recirculation  zone  length  correlation  used  here  was  based  on 
Drewry's  results  [l 6 ] .  Previous  work  also  showed  that  the  dividing  stream¬ 
line  shear  stress  coefficient  had  a  strong  effect  on  the  predicted  static 
pressure  distribution,  while  the  turbulent  kinetic  energy  dissipation  rate 
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length  scale  assumed  for  the  recirculation  zone  had  a  strong  effect  on 
the  rate  of  transport  into  and  out  of  the  recirculation  region.  The 
shear  layer  initial  thickness  and  rate  of  growth  had  little  overall 
effect,  over  the  range  considered  in  the  earlier  work. 

Although  Chaturvedi's  data  show  a  reattachment  point  at  9.2 
step  heights  from  the  inlet  and  indicate  an  essentially  parabolic  re¬ 
circulation  zone  shape,  comparison  of  the  predictions  of  the  model  with 
the  data  for  static  pressure  and  centerline  velocity  showed  that  these 
data  could  not  be  reproduced  with  the  parabolic  recirculation  zone 
boundary  that  has  been  used  in  other  modular  model  calculations.  In¬ 
stead,  an  extended  recirculation  zone  shape  had  to  be  adopted.  Figure  5 
illustrates  the  experimental ly  determined  reci rculation  zone  shape,  and 
the  parabolic  and  modified  recirculation  zone  shapes  adopted  for  the 
computation.  Note  that  when  the  extended  shape  is  used,  the  effective 
recirculation  zone  length  for  purposes  of  computing  the  stirred  reactor 
volume  is  kept  at  the  1.61  ft.  distance  indicated  by  the  data;  thus  the 
extension  represents  a  crude  model  of  a  region  over  which  the  displace¬ 
ment  effect  of  the  shear  layer  bounding  the  dividing  streamline  relaxes 
to  that  of  a  wall  boundary  layer. 


A  DATA,  CHATURVEDI [l5] 

-  PARABOLIC  DIVIDING  STREAMLINE  MODEL 

-  MODIFIED  DIVIDING  STREAMLINE  MODEL 


1.0  2.0 

AXIAL  DISTANCE,  FT. 
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Figure  5.  COMPARISON  OF  ASSUMED  DIVIDING  STREAMLINE  SHAPES 
WITH  CHATURVEDI  DATA  [ 15] 
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Centerline  velocity  and  turbulent  kinetic  energy  profiles  and 
wall  static  pressure  distributions  are  shown  in  Figure  6  in  comparison 
with  Chaturvedi's  data  [15].  For  these  calculations,  the  initial  velocity 
was  taken  to  be  uniform  across  the  inlet,  at  100  ft/sec,  the  initial  tur¬ 
bulent  kinetic  energy  level  was  10_4Uq2  ,  based  on  Chaturvedi's  data,  and 
the  initial  turbulence  energy  dissipation  rate  was  obtained  from  e  = 

1.69k1  •  s/i-|c  with  5,^  =  0.1  rg.  A  standard  set  of  two-equation  model  coef¬ 
ficients  have  been  used  for  all  of  the  flowfield  computations  reported  in 
this  paper;  thus: 

C^  =  0.09;  CE1  =  1.40;  C£2  =  1.95;  oR  =  1.00;  =  1.22. 

Results  are  shown  for  two  values  of  the  effective  shear  stress 
coefficient  along  the  dividing  streamline,  with  CFW  "  0  .016  providing  the 
best  agreement  with  the  data  for  wall  static  pressure  and  centerline 
velocity,  but  underpredicting  the  centerline  kinetic  energy  data.  Note 
that  the  initial  increase  of  the  centerline  kinetic  energy  is  not  predic¬ 
ted  by  the  turbulence  model  used  here;  this  phenomenon,  which  is  associa¬ 
ted  with  the  closure  of  the  potential  core  region  of  the  flow,  has  also 
been  observed  in  free  jets.  The  other  model  parameters  used  in  this  cal¬ 
culation  are  a  =  0.03,  b  =  0.09,  kR2  =  0.005  Ug2,  and  £RZ  =  (h  +  LR?)/2, 
where  the  latter  two  variables  represent  the  effective  turbulent  kinetic 
energy  and  dissipation  length  scale  within  the  recirculation  region.  In 
the  expression  for  the  recirculation  zone  length  scale,  h  is  the  step 
height  and  LRZ  the  recirculation  zone  length;  since  LR^  =  9.2  h,  effec¬ 
tively  £R£  =  5.1  h. 

The  predicted  velocity  profiles  are  compared  with  the  experi¬ 
mental  data  on  Figures  7  and  8.  Overall,  the  agreement  is  reasonably 
good,  although  it  is  clear  that  the  flowfield  mixing  rate  is  slightly 
overpredicted.  For  X/D  =  2  and  X/D  =  4,  the  predicted  velocity  profiles 
end  at  the  dividing  streamline  location,  since  the  modular  model  does 
not  provide  any  detail  of  the  recirculation  region  flowfield.  Figure  9 
shows  a  comparison  of  predicted  radial  turbulent  kinetic  energy  profiles 
at  axial  locations  downstream  of  the  recirculation  zone;  again  the  overall 
level  of  agreement  is  reasonably  good. 
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This  comparison  provides  some  validation  of  the  ability  of  the 
modular  approach  to  model  accurately  the  details  of  a  given  flowfield. 

It  is  important  to  note  that  the  adjustment  of  the  model  parameters  car¬ 
ried  out  for  this  case  was  done  with  reference  to  the  wall  static  pres¬ 
sure  and  centerline  velocity  measurements :  the  detailed  profile  data 
were  not  used  in  the  modeling  iteration.  In  the  course  of  the  modeling 
of  this  flowfield,  a  number  of  observations  of  model  sensitivities  were 
made: 

•  the  predictions  wore  insensitive  to  the  parameters  a  and 
b  in  the  shear  layer  growth  model; 

•  the  initial  level  of  turbulent  kinetic  energy  had  a  strong 
effect  on  the  overall  turbulent  kinetic  energy  distribution 
in  the  flowfield,  but  a  significantly  smaller  effect  on  the 
wall  pressure  distribution; 

•  the  wall  pressure  distribution  is  primarily  influenced  by 
the  size  and  shape  of  the  recirculation  zone  and  by  the 
level  of  the  assumed  dividing  streamline  shear  stress 
coefficient,  CpW. 

It  was  noted  earlier  that  the  incompressible  sudden-expansion 
represents  a  conf iguration  for  which  the  computational  capabilities  of 
the  modular  and  unified  models  overlap,  so  that  for  this  flowfield  a 
comparison  of  results  obtained  using  both  types  of  model  is  of  consider¬ 
able  interest.  To  provide  this  comparison,  calculations  of  the  Chaturvedi 
sudden  expansion  flowfield  were  carried  out  using  an  elliptic  code  de¬ 
veloped  from  the  Imperial  College  primitive  variables  formulation  [17]. 

The  inlet  conditions  and  turbulence  modeling  used  for  the  elliptic  cal¬ 
culation  were  identical  to  that  used  in  the  modular  approach;  wall  bound¬ 
ary  conditions  for  the  elliptic  computations  were  defined  using  a  "wall- 
function"  technique  [ 1 7]  and  the  downstream  boundary  condition  was  ob¬ 
tained  through  use  of  a  zero  axial  gradient  hypothesis  for  all  dependent 
variables.  The  unified  model  calculation  used  a  100  x  25  (x  -  r)  grid 
to  minimize  numerical  diffusion  problems,  and  was  continued  for  1200 
iterations  to  obtain  a  "mass  imbalance"  convergence  criterion  of  less 
than  0.01%. 

Comparisons  of  radial  velocity  profiles  obtained  through  use 
of  the  unified  model  with  those  obtained  from  the  modular  calculation 


and  with  Chaturvodi's  [lo]  experimental  data  are  shown  in  Figures  10  and 
11.  It  is  evident  from  these  figures  that  the  unified  model  predicts  a 
greater  mixing  rate  and  shorter  recirculation  zone  length  than  the  data 
indicate;  however,  no  effort  was  made  to  optimize  the  unified  model  cal¬ 
culations  for  these  data  (through,  for  example,  modifications  to  the  wall 
boundary  conditions  used  in  the  computation).  It  is  interesting  to  note 
that  for  X/D  =  2  and  X/D  =  4,  the  unified  model  velocity  profile  shape 
predictions  are  similar  to  those  the  modular  model  produces,  up  to  the 
dividing  streamline.  The  modular  model ,  of  course,  provides  no  detail 
within  the  recirculation  region.  Also  of  interest  is  the  steepening  of 
the  velocity  profile  in  the  near  wall  region  between  X/D  =  6  and  X/D  =  8 
that  is  predicted  by  both  models.  Figure  11.  This  provides  further  evi¬ 
dence  that  the  shear  layer  relaxation  region  hypothesized  to  exist  in  the 
modular  model  formulation  of  the  flowfield  is  in  fact  a  feature  of  this 
flow. 
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Figure  10.  COMPARISON  OF  VELOCITY  PROFILES  PREDICTED  BY 
MODULAR  AND  UNIFIED  MODELS  WITH  DATA  OF 
CHATURVEDI  [15] 


Figure  11.  COMPARISON  OF  VELOCITY  PROFILES  PREDICTED  BY 
MODULAR  AND  UNIFIED  MODELS  WITH  DATA  OF 
CHATURVF.DI  [15] 

Radial  turbulent  kinetic  energy  profiles  for  X/D  =  6  are  shown 
in  Figure  12.  Although  the  elliptic  model  produces  a  slightly  lower 
level  of  turbulent  kinetic  energy  than  does  the  modular  approach,  both 
calculations  are  in  reasonably  good  agreement  with  the  data. 

These  comparisons  indicate  that  for  this  configuration,  the 
modular  approach  is  capable  of  providing  predictions  of  the  details  of 
the  flowfieid  that  are  in  agreement  with  both  the  experimental  data  and 
with  predictions  obtained  from  unified  techniques.  Although,  as  has  been 
noted  previously,  the  modular  model  is  not  a  predictive  technique  in  the 
sense  that  unified  approaches  are  intended  to  be,  these  comparisons  indi¬ 
cate  that  the  shear  layer  and  recirculation  zone  modeling  inherent  in  the 
modular  approach  are  both  adequate  representations  of  these  regions. 

This  provides  further  validation  of  the  modular  model  approach  which,  for 
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reacting  sudden  expansion  flows,  is  capable  of  providing  far  more  detail 
than  is  practically  available  with  unified  techniques. 


r 

rl 

Figure  12.  COMPARISON  OF  TURBULENT  KINETIC  ENERGY  PROFILES 
PREDICTED  BY  MODULAR  AND  UNIFIED  MODELS  WITH  DATA 
OF  CHATURVEDI  [15] 

The  configuration  used  experimentally  by  Craig,  et  al .  [18]  is, 
while  still  a  sudden-expansion,  considerably  different  from  that  studied 
by  Chaturvedi  [l5j.  The  basic  configuration  studied  in  reference  18  in¬ 
volved  a  combustor  L/D  of  3,  a  dump  area  ratio  A3/A2  of  2.25,  and  an  exit 
nozzle  area  ratio  A^/A^  -  0.40.  Thus  it  is  considerably  shorter  than 
Chaturvedi 's  apparatus,  which  had  an  L/D  of  13  based  on  the  "combustor" 
diameter,  and  incorporates  an  exit  nozzle  contraction  which  was  not  used 
in  Chaturvedi' s  work.  Doth  of  those  differences  can  be  expected  to  have 
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substantial  effects  on  the  overall  flowfield  development.  Two  sets  of 
inlet  ain  low  rates  and  inlet  total  temperatures  w-ro  used  throughout  the 
wort  reported  in  reference  18,  with  nominal  values  of  1.4b  kg/sec  and 
0.73  kg/sec,  and  bbb°K  and  722,0K,  respectively.  A  major  contribution  of 
this  work  is  the  measurement  of  wall  static  pressure  distributions  in 
both  cold  flow  and  reacting  flow.  This  report  represents  one  of  the  few 
investigations  in  which  these  data,  crucial  to  the  comparison  of  experi¬ 
mental  data  with  analytical  model  results,  have  been  reported.  Supple¬ 
menting  tie  results  reported  in  reference  18,  additional  unreported 
static  pressure  distributions  have  been  made  available  to  the  authors  by 
Dr.  Craig. 

Although  an  appropriate  set  of  model  coefficients  have  been 
develop'd  for  Chaturvedi’s  configuration,  the  differences  between  that 
apparatus  and  the  combustor  used  by  Craig,  et  al .  are  substantial.  Thus 
it  can  Ir-  expected  that  different  parameter  values  would  be  necessary  to 
accurately  characterize  Craig's  combustor  flowfield.  In  order  to  estab¬ 
lish  the  necessary  values,  parametric  calculations  of  one  set  of  cold 
flow  data,  for  an  inlet  total  temperature  of  552°K  and  an  air  mass  flow 
of  1.45  l.f t/scc ,  were  carried  out.  Figure  13  illustrates  the  results  of 
these  co:  putations.  Here,  a  parabolic  recirculation  zone  shape  was  used, 
and  those  model  parameters  not  noted  on  Figure  13  retained  the  values 
that  wore  arrived  at  in  the  course  of  the  computations  of  Chaturvedi's 
configuration.  The  inlet  turbulent  kinetic  energy  level  and  dissipation 
rate  were  established  using  the  same  relationships  to  mean  flow  and  geo¬ 
metric  variables  as  were  used  in  the  Chaturvedi  calculations. 

The  differences  in  model  parameters  between  the  Chaturvedi  [lb] 
and  Craig,  et  al.  [l8]  cold  flow  configurations  are  the  shape  of  the  re¬ 
circulation  zone  and  the  level  of  the  dividing  streamline  shear  stress 
coefficient.  Although  the  comparisons  shown  in  Figure  13  may  indicate 
that  a  slight  further  refinement  to  the  parabolic  shape  used  in  these 
calculations  could  improve  the  match  between  the  computed  and  measured 
pressures  over  the  range  0.5  <  X/D  <  1.5,  overall  the  level  of  agreement 
is  quite  good.  The  apparent  lack  of  a  large  transition  region  between 
shear  layer  and  wall  boundary  layer  implied  by  these  results,  obtained 
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with  no  tr  nsi  tion  region,  may  be  u  resul  t.  of  'ho  use  of  a  convergin  j 
exit  nn;'/i ,  ■  in  the  Crain,  et  al .  confirm ti<  o.  A.  fuo'nor  difference 
between  the  Craig,  et  al  .  configuration  uni  that  studf".!  try  Cliat urvedi 
is  the  lath  of  an  inlet  section  upstream  of  ».he  tiir.p  pi. me  in  the  latter 
case,  coa.j  ired  to  the  use  of  a  relatively  long  inlet  section  in  *  he 
for  er.  Tiiese  computations  use  a  uniform  inlet  profile  so  that  the 
existence  of  an  inlet  boundary  layer  is  reflected  primarily  in  the  shear 
stress  pvu-eter  Cp^.  Thus  the  requirement  that  tne  skin  friction  coef¬ 
ficient  along  the  dividing  streamline  be  a  factor  of  3  higher  for  the 
Craig,  c?  al.  comparison  than  for  the  Chaturvedi  modeling  may  be  a  re¬ 
sult  of  the  differences  in  the  initial  conditions  between  the  two  flows. 
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Figure  13.  COMPARISON  OF  PREDICTED  WALL  STATIC 

PRESSURE  DISTRIBUTIONS  WITH  COLD  FLOW 
SUDDEN -EXPANSION  COMBUSTOR  DATA  [ 18] 

3.2  REACTING  FLOW  APPLICATIONS 

Craig,  et  al .  [18]  also  report  combustion  efficiency  and  wall 
static  pressure  distribution  measurements  for  a  premixed  sudden  expansion 
combustor.  The  combustor  configuration  is  identical  to  that,  used  in  the 
cold  flow  investigation.  To  further  test  the  modular  model  formulation, 
computations  were  made  of  a  premixed  dump  combustor  at  a  fuel -air  ratio 
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of  0.0f>3  for  comparison  wifi'  the  data  presented  in  reference  18.  In  the 
experiments,  the  fuel  used  v.ms  J!'-4  which  was  represented  in  the  computa¬ 
tions  by  propane.  Both  a  simple  one-step  global  finite  rate  chemistry 
model  and  the  full  hydrocarbon  oxidation  kinetics  as  represented  by  the 
quasiglobal  model  [c,  19]  were  used  to  represent  the  chemical  kinetics 
processes  in  the  directed  flow.  The  global  model  was  used  for  the  recir¬ 
culation  zone  in  both  cases.  Propane  was  chosen  to  represent  the  fuel 
because  previous  studies  with  propane  had  resulted  in  the  development  of 
a  one-step  finite  rate  kinetics  model  that  represented  fairly  accurately 
the  ignition  delay  for  propane-air  (although  not  the  overall  reaction 
time)  over  a  range  of  conditions  of  interest  in  this  work. 

Initial  conditions  included  an  inlet  total  temperature  of  554°K 
(997°R)  and  an  inlet  static  pressure  of  1.83  atm  (3881  lb/ft2).  With  a 
mass  flow  rate  of  1.87  kg/sec  (3.49  lb/sec),  the  inlet  velocity  and  static 
temperature  were  199  m/sec  (931  ft/sec)  and  543°K  (977°R),  respecti vely; 
the  inlet  Mach  number  was  0.351.  Initial  turbulent  kinetic  energy  and 
dissipation  rate  values  were  established  in  the  same  manner  as  for  the 
cold  flow  calculations ,  and  the  geometry  of  the  combustor  and  recircula¬ 
tion  zone  was  the  same  as  in  the  cold  flow  calculations. 

Both  combustion  efficiency  and  wall  static  pressure  distribu¬ 
tion  data  arc  available  for  this  configuration,  and  the  results  of  the 
modular  model  calculation  of  those  quantities  are  shown  in  comparison 
with  the  experimental  data  in  Figure  14.  The  combustion  efficiency  shown 
was  computed  from  the  mass-average  temperature  at  each  axial  location  in 
the  calculation,  using  the  JANNAF  temperature-rise  combustion  efficiency 
definition,  reference  ?0*.  The  "ideal"  exit  total  temperature  required 
for  this  calculation  was  obtained  from  the  data  tabulations  provided  to 
the  authors  by  Hr.  R.  R.  Craig  of  AF APL ;  for  these  data,  Tj^(ideal)  - 
2178°K.  Note  that  in  these  comparisons,  data  points  are  shown  for  the 
fuel/air  ratio  0.093  value  at  which  the  computations  were  carried  out. 
These  data  points  were  obtained  from  plotted  values  of  combustion  effi¬ 
ciency  vs.  fuel -air  ratio  presented  in  reference  18  for  three  different 

*  The  de'inition  was  modified  slightly  in  that,  static  rather  than  total 
ten-no ra t ures  were  used,  for  convenience.  At  the  combustor  Mach  numbers 
emoentiisd  is  Mi«’se  experiim  nts ,  the  difference  in  temperatures  is  of 
the  ord -  i-  of  i  .  , 
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values  or  r.o:;ibustor  L/0.  The  band  shown  for  each  data  point  represents 
the  ratine  of  observed  combustion  efficiencies  as.  a  function  of  fuel-air 
ratio  and  is  a  better  indication  of  the  overall  trend  of  the  combustion 
efficiency  vs.  fuel-air  ratio  data  than  are  the  individual  data  points 
themselves.  Some  caution  is  advised  in  interpreting  the  combustion  effi¬ 
ciency  comparison  shown  in  figure  14  since  three  different  combustor  con¬ 
figuration.  were  involved  in  obtaining  the  data.  Thus  the  relationship 
between  recirculation  zone  length  and  combustor  length  is  different  for 
each  of  the  three  combustors  tested.  On  the  other  hand,  the  static  pres¬ 
sure  data  (for  two  fuel -air  ratios  which  bracket  the  fuel -air  ratio  used 
in  the  computation)  and  the  predictions  are  both  for  a  combustor  L/D  of 
3.0. 

The  first  point  that  should  be  noted  in  reviewing  the  results 
presented  in  Figure  14  is  that  the  level  of  agreement  with  the  data  ob¬ 
tained  was  arrived  at  by  increasing  the  kinetic  energy  level  in  the  re¬ 
circulation  zone  substantially  over  the  levels  found  to  be  appropriate 
for  the  cold  flow  cases  already  discussed.  One  possible  physical  cause 
for  this  increased  turbulent  energy  level  in  the  recirculation  region  is 
the  large-scale  oscillation  of  the  recirculation  zone  that  has  been  ob¬ 
served  in  a  variety  of  reacting,  recirculating  flowfields  [19].  While 
the  fluctuations  that  such  oscillations  produce  are  not  strictly  turbu¬ 
lence,  two  factors  combine  to  suggest  that  they  can  be  interpreted  in  the 
context  of  an  analytical  flow  model  as  increased  turbulent  intensity 
levels  in  the  recirculation  zone.  One  is  that  the  oscillations  enhance 
the  rate  of  mixing  in  the  recirculation  region  and  in  the  shear  layer 
bounding  it,  and  the  second  is  that  some  of  the  energy  in  these  oscilla¬ 
tions  can  be  expected  to  provide  a  source  of  turbulence  energy.  For  all 
other  fartors  equal,  the  increase  in  turbulent  kinetic  energy  in  the 
model  recirculation  zone,  while  not  increasing  the  mixing  rate  within 
that  region,  does  increase  the  effective  viscosity  along  the  dividing 
streamline,  and  thus  the  rate  of  transport  across  the  shear  layer.  It 
also  changes,  both  directly  and  indirectly  through  the  transport  mecha¬ 
nism,  the  rate  of  transport  of  turbulent  kinetic  energy  to  the  recircula¬ 
tion  zone.  Thus  the  increased  turbulent  energy  level  in  the  recirculation 
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zone  required  to  achieve  cjood  agreement  with  the  available  data  for  the 
Craig,  et  al.  premixed  combustor  appears  to  be  plausible  on  physical 
grounds . 
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Figure  14.  COMPARISON  Of  PREDICTED  AND  MEASURED 
COMBUSTION  EFFICIENCY  AND  WAIL  STATIC 
PRESSURE  DISTRIBUTIONS,  PREMIXED  SUDDEN- 
EXPANSION  COMBUSTOR  GLOBAL  FINITE  RATE 
CHEMISTRY  MODEL 
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It  is  I  so  of  interest  to  note,  from  Figure*  14,  that  1 1 1 « •  lower 
of  the  two  curves  shown  in  the  combustion  efficiency  comparison  provides 
the  better  agreement  with  the  observed  static  pressure  profiles,  despite 
the  significant  underprediction  by  this  computation  of  the  combustion 
efficiency  at  X/D  -  3.0.  As  has  already  been  noted,  the  global  finite 
rate  model  does  not  properly  predict  reaction  time  for  the  propane-air 
system  used  to  model  this  flow,  and  thus  further  computations  were  car¬ 
ried  out  using  the  quasiglohal  hydrocarbon  kinetics  model  [C,  19]  for 
the  same  conditions. 

The  initial  quasiglohal  computation  of  this  flowfield,  denoted 
by  F  -  1  on  Figure  15,  produced  a  surprising  result:  blowout  was  predic¬ 
ted  to  occur  in  the  directed  flow.  To  interpret  this  result,  a  digression 
to  a  discussion  of  the  development  of  the  quasiglohal  chemical  kinetics 
model  for  hydrocarbon  oxidation  is  required. 

During  the  development  of  the  quasiglohal  model,  two  major 
sources  of  data  were  utilized.  One  of  these  was  ignition  delay  data  ob¬ 
tained  in  experimental  configurations  approximating  the  ideal  limit  of 
plug  flow,  and  the  other  was  blowout  data  obtained  from  laboratory  well- 
stirred  reactors.  Although  data  from  both  sources  verified  the  basic 
quasiglohal  model  concept:  that  hydrocarbons  higher  than  propane  could 
be  grouped  into  generic  classifications  and  that  for  each  group  an  over¬ 
all  partial  oxidation  reaction  to  the  products  CO  and  could  be  written, 
at  least  for  fuel-air  equivalence  ratios  less  than  one,  it  was  observed 
that  the  effective  pre-exponential  coefficient  in  the  quasiglobal  Ahrren- 
ius  rate  expression  had  to  be  increased  by  a  factor  of  80  for  well-stirred 
reactor  flows  compared  to  plug  flow. 

The  salient  difference  between  a  well-stirred  reactor  and  a  plug 
flow  is  that  in  the  former  the  reactants  and  products  of  combustion  are 
continuously  being  mixed  with  fresh,  unreacted,  entering  species.  This, 
of  course,  is  the  situation  for  a  mixing  and  reacting  flow,  and  it  is  ex¬ 
pected  to  he  particularly  significant  near  blowout.  Thus,  a  further 
series  of  computations  were  carried  out  with  a  multiplicative  factor  F 
preceding  the  pre-exponential  coefficient  in  the  quasiglobal  hydrocarbon 
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part ial  oxidation  rate  expression  [19].  For  propane,  this  expression  is 

kf  -  6.0  x  10“  T  exp  (-12,?8b/T)  (42) 

and  values  of  F  from  the  plug  flow  limit  of  1  to  the  stirred  reactor 
limit  of  tit)  were  considered.  The  rationale  for  carrying  out  this  investi¬ 
gation  is  that  each  element  of  reacting  fluid  can  be  considered  to  behave 
as  a  combination  of  plug  flow  and  stirred  reactor  flow:  in  the  absence 
of  turbulent  transport  each  element  is  a  plug  flow  reactor,  whereas  if 
turbulent  transport  dominates  then  the  flow  in  each  element  approaches 
that  of  a  homogeneous  stirred  reactor.  This  idealization  is  of  course  an 
oversimp] ification  but  it  is  an  instructive  concept. 

As  is  shown  in  Figure  15,  the  well-stirred  reactor  limit,  F  = 

80,  drastically  overpredicts  the  rate  of  reaction  in  comparison  with  the 
data,  while  for  F  =  5,  blowout  is  still  predicted.  Wh.cn  F  =  10  is  selec¬ 
ted,  an  almost  perfect  representation  of  the  data  is  obtained,  but  this 
should  In  considered  to  be  fortuitous:  no  representation  that  F  =  10  is 
the  appiv  riate  value  for  mixing  flows  can  or  should  be  made. 

If  the  predicted  pressure  profiles  are  considered,  figure  15, 
an  entirely  different  conclusion  could  be  obtained.  Examination  of  these 
predictions  shows  that,  while  none  of  the  predictions  is  in  exact  corre¬ 
spondence  with  the  data,  the  predictions  for  F  -  1  and  F  =  5  are  in  better 
agreement  than  those  for  F  =  10  and  F  =  80.  Significantly,  the  disagree¬ 
ment  in  the  latter  two  cases  is  greatest  in  the  region  downstream  of  the 
end  of  the  reci rculation  zone  for  this  configuration,  and  is  thus  not 
related  to  the  dividing  streamline  modeling  incorporated  in  the  model 
approach.  The  data  presented  by  Craig,  et  al .  [l8]  indicate  that  for  the 
premixed  combustor,  blowout  is  observed  between  a  fuel-air  ratio  of  0.030 
and  0,040.  The  sensitivity  of  these  results  to  the  parameter  F,  observed 
in  Figure  15,  is  an  indication  that  at  the  conditions  of  these  data,  a 
fuel-air  ratio  jf  0.053,  the  flowfield  is  also  near  blowout.  In  fact, 
the  combination  of  the  agreement  with  the  axial  pressure  profiles  obtained 
for  F  1  and  F  -  5  and  the  agreement  with  the  measured  combustion  effi¬ 
ciency  for  F  ~  10  can  he  taken  to  indicate  that  this  combustor  flow  was, 
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in  fart,  oscillatory,  with  a  continuous  partial  blowout  and  reignition 
phenomena  occurring  during  the  test. 
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Fiqure  15.  COMPARISON  OF  PREDICTED  AND  MEASURED 
COMBUSTION  EFFICIENCY  AND  WALL  STATIC 
PRESSURE  DISTRIBUTIONS;  DATA  FROM 
CRAIG,  et  al .  [18] 


This  discussion  has  emphasized  the  model  sensitivities  observed 
during  this  phase  of  the  model  validation  program.  It  is,  of  course, 
also  instructive  to  consider  the  insensitivities  observed.  These  can  be 
summarized  as  follows: 

•  For  the  premixed  reacting  flow,  the  predictions  are  insensi¬ 
tive  to  the  details  of  the  shear  layer  growth  rate,  as 
modeled  by  Eq.  29. 

•  The  premixed  combustor  is  totally  insensitive  to  initial 
velocity  profile  variation  as  represented  by  the  inclu¬ 
sion  in  the  initial  profile  of  a  1/7  power  law  boundary 
layer  with  a  thickness  equal  to  0.1  times  the  inlet  radius. 
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SECTION  4 

APPLICATION  OF  UNIT  ANALYSES  TO  COMBUSTOR  PROBLEMS 


Although  this  report  has  focused  on  the  application  of  the 
complete  modular  model  to  the  analysis  and  interpretation  of  combustor 
data,  unit  analyses  can  also  provide  considerable  insight  into  various 
aspects  of  combustor  performance.  These  unit  analyses  can  themselves 
be  elements  of  an  overall  modular  combustor  model:  one  of  the  advantages 
of  modular  modeling  is  that  the  model  elements  can  be  individually  devel¬ 
oped  and  verified  through  use  as  unit  analysis  tools.  In  this  section, 
two  types  of  unit  analyses  will  be  discussed.  The  first  example  involves 
the  use  of  a  simple  droplet  vaporization  analysis  to  interpret  experimen¬ 
tal  data  in  a  laboratory  combustor,  while  the  second  example  involves  the 
use  of  the  well -stirred  reactor  module  of  the  modular  model  as  a  unit 
analysis  of  flameholder  blowout  phenomena. 

4.1  DROPLET  VAPORIZATION  EFFECTS  ON  COMBUSTOR  PERFORMANCE 

Because  of  the  requirement  for  greater  energy  density  than 
conventional  fuels,  ramjet  fuels  can  have  considerably  different  physical 
characteristics  from  fuels  such  as  OP-4.  These  physical  characteristics 
include  viscosity,  which  is  generally  greater  for  ramjet  fuels,  and  vola¬ 
tility,  which  is  in  general  lower.  These  differences  can  be  expected  to 
lead  to  larger  and  more  persistent  fuel  droplets  for  ramjet  fuels  than 
for  conventional  fuels. 

The  combustion  characteristics  of  several  ramjet  fuels  were 
investigated  by  Siminski  [ 2 1 ] ,  who  presented  data  for  the  specific  im¬ 
pulse  (ISP)  efficiency  obtained  using  various  fuels  in  a  simulated  ramjet 
combustor.  These  data  show  a  strong  dependence  of  n j on  inlet  tempera¬ 
ture  for  low  inlot  temperatures.  The  fuel  was  injected  into  the  combus¬ 
tor  in  liquid  form,  so  that  the  vaporization  rate  of  the  fuel  can  be 
expected  to  have  an  impact  on  the  overall  combustor  performance.  Since 
njt-p  is  related  to  the  actual  total  temperature  rise  in  the  combustion 
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process  divided  by  the  ideal  total  temperature  rise,  and  this  can  in 
turn  be  related  to  the  vapor  fuel  concentration  divided  by  the  initial 
liquid  fuel  concentration,  which  in  turn  is  directly  relatable  to  the 
fuel  vaporization  rate,  a  relation  between  n^p  and  vaporization  rate 
can  be  devised: 


and 


‘ISP 


Tq  actual \2 
Tg  ideal  ) 


T  ,0  / 


W; 

(43) 

This  correlation  was  tested  by  computing  vaporization  rates  as  a  function 
of  gas-droplet  AT  for  two  droplet  sizes  and  for  two  fuels,  RJ-5  and  JP-6, 
which  differ  substantial ly  in  viscosity  and  boiling  point  temperature.  A 
simple  one-dimensional  droplet  vaporization  computation  was  used,  incor¬ 
porating  the  "d2"  vaporization  law 


d2 


(44) 


where  Nu  is  the  Nusselt  number  and  X  is  an  evaporation  rate  constant 
which  depends  on  droplet  thermal  conductivity,  density,  specific  heat 
and  heat  of  vaporization  and  on  the  droplet-gas  phase  temperature  differ¬ 
ence.  The  solution  for  droplet  diameter,  d,  is  obtained  through  use  of 
the  "d2"  law  coupled  to  one-dimensional  statements  of  momentum  and  energy 
transport  for  the  droplet  exposed  to  a  uniform  gas-phase  environment. 
Droplet  diameter  was  computed  at  a  time  from  injection  corresponding 
approximately  to  the  time  required  for  the  fuel  to  reach  the  flameholder 
in  the  combustor  used  by  Siminski  [Zl].  The  agreement  between  the  corre¬ 
lation  and  the  data  is  striking  for  both  fuels,  as  shown  in  Figures  16 
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and  17.  Since  the  droplet  size  in  the  apparatus  used  by  S i mi  ns  k i  was 
not  known,  a  range  of  sizes  was  used;  it  might  be  expected  that  at 
similar  conditions  the  RJ-5  would  form  larger  droplets  than  the  JP-5 
because  of  its  higher  viscosity.  Note  also  from  figure  16  that  the 
agreement  of  the  correlation  for  smaller  droplets  with  the  data  at  the 
higher  fuel  temperature  is  consistent  with  a  reduction  in  drop  size  with 
fuel  temperature  increase  (due  to  a  fuel  viscosity  increase);  this  is 
also  seen  for  JP-5,  Figure  17  for  the  lower  AT  (and  therefore  higher  fuel 
temperature)  range. 


nI$P 

Figure  16.  CORRELATION  BETWEEN  VAPORIZATION  RATE  AND 
ISP  EFFICIENCY,  RJ-5  FULL 
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Figure  17.  CORRELATION  BETWEEN  VAPORIZATION  RATE 
AND  ISP  EFFICIENCY,  JP-5  FUEL 


4.2  FLAMEHOLDER  BLOWOUT  AT  HIGH  INLET  TEMPERATURE 

The  well-stirred  reactor  analysis  is  a  model  of  a  laboratory 
well-stirred  reactor.  It  is  also  a  model  of  a  recirculation  region, 
such  as  that  downstream  of  the  dump  station  in  a  sudden-expansion  com¬ 
bustor,  or  that  behind  a  flaneholder.  Thus  the  well-stirred  reactor 
model  can  be  used  as  a  unit  analysis  of  flameholder  flame  stabilisation 
phenomena . 
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An  example  of  the  use  of  the  well-stirred  reaclor  model  in  the 
study  of  flame  stabilization  phenomena  is  the  study  of  the  effects  of 
vitiation  on  flameholder  performance.  In  order  to  achieve  high  combustor 
inlet  temperatures  in  a  ramjet  ground  test  facility,  vitiation  of  the 
inlet  air  stream  is  often  necessary,  and  even  if  a  hydrocarbon  fuel  such 
as  propane  is  used  as  the  vitiator  fuel,  with  oxygen  makeup,  serving  to 
minimize  potential  molecular  weight  mismatch  between  the  vitiated  test 
stream  and  the  inlet  air  stream  at  flight  conditions,  the  question  of  the 
effects  of  vitiation  on  measured  flameholder  performance  still  arises. 
Flame  stabilization  can  be  modeled  using  the  well-stirred  reactor  as  an 
approximation  to  an  intensely  backmixed  recirculation  region  behind  a 
flameholder.  The  loading  parameter  for  this  well-stirred  reactor  at 
blowout  then  becomes  a  characterization  of  the  flame  stability  which  can 
be  expected.  To  study  the  effects  of  vitiation  on  flame  stabilization 
at  high  inlet  temperatures,  well-stirred  reactor  calculations  were  car¬ 
ried  out,  using  both  pure  ana  vitiated  air  inlet  streams. 

Results  of  these  computations  are  shown  in  Figure  18  for  one 
test  condition  of  interest.  The  effect  of  vitiation  for  these  conditions 
can  be  seen  to  be  a  reduction  of  the  loading  parameter  at  blowoff,  because 
of  the  effects  of  the  presence  of  H^O  and  CO^  on  the  specific  heat  of  the 
"air"  mixture. 

These  results  are  presented  in  terms  of  a  reduced  loading  param¬ 
eter  as  used  in  Curran  [22];  in  general  the  reduced  loading  parameter  is 
defined  as 


(45) 


with  A  and  b  designated  as  empirical  parameters.  In  actuality,  the  expo¬ 
nent  on  the  pressure  is  itself  an  empirical  parameter.  If  the  non-vitia- 
ted  results  shown  in  Figure  18  are  considered,  a  comparison  of  these  data 
with  the  available  experimental  results  (from  references  22,  23  and  19) 
indicates  that,  for  a  pressure  exponent  of  1.8,  the  temperature  exponent 
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b  varies  from  1.0  at  an  equivalence  ratio  of  0.4  to  1.46  at  an  equivalence 
ratio  of  1.0,  as  shown  by  the  additional  points  on  figure  18.  However, 
the  data  shown  in  Figure  18  were  obtained  at  one  atmosphere  pressure, 
while  the  computations  are  for  a  pressure  level  of  about  two  atmospheres. 
Thus  the  effects  of  the  use  of  a  pressure  exponent  of  1.8  are  also  in¬ 
cluded  in  these  results.  Further  computations,  at  temperature  and  pres¬ 
sure  levels  that  are  elevated  with  respect  to  the  available  data,  can 
serve  to  indicate  the  appropriate  values  for  both  the  pressure  and  tem¬ 
perature  exponents  in  Eq.  45. 
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Figure  18.  COMPARISON  OF  COMPUTED  BLOWOUT  RESULTS  FOR 

HIGH  TEMPERATURE  CONDITIONS  WITH  EXPERIMENTAL 
STIRRED  REACTOR  DATA 
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SECTION  5  -  CONCLUSIONS 


While  the  computational  ability  to  predict  in  detail  the 
flowfield  in  a  ramjet  combustor  without  empirical  approximations  appears 
to  remain  a  long-term  goal  rather  than  a  short-term  prospect,  modeling 
techniques  have  been  shown  to  be  a  viable  tool  in  the  analysis  of  ramjet 
combustion  processes.  These  modeling  techniques  range  from  the  analysis 
of  plug  flow  and  well-stirred  reactor  flows,  which  represent  idealiza¬ 
tions  of  aspects  of  an  overall  ramjet  flowfield,  to  the  modular  approach, 
which  represents  the  only  realistic  method  of  introducing  a  detailed 
description  of  a  ramjet  combustor  flowfield  into  a  computational  technique. 

The  modular  approach  is  not  a  predictive  technique,  and  no  cur¬ 
rently  available  mathematical  model  of  a  ramjet  combustor  can  be  considered 
to  be  truly  predictive.  What  the  modular  approach,  coupled  with  the  use 
of  detailed  unit  analyses  offers,  is  the  ability  to  analyze  and  interpret 
experimental  data  and  to  aid  in  the  design  and  planning  of  appropriate 
experiments . 

A  considerable  amount  of  further  work  is  still  required  to 
develop  all  viable  ramjet  modeling  approaches,  whether  unified  or  modular 
formulations.  What  the  work  described  in  this  report  has  shown  is  the 
insight  that  can  be  gained  through  t.hc  use  of  modeling  in  both  the  design 
of  combustors  and  experiments  and  in  the  interpretation  of  experimental 
results.  It  is  in  this  area  that  modeling  holds  the  greatest  promise. 
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APPENDIX  I 

NUMERICAL  SOLUTION  TECHNIQUE  FGR  DIRECTED  PLOW  ANALYSIS 


The  solution  technique  used  Tor  the  directed  Flow  portion  of  the 
modular  model  involves  the  use  of  a  Von-Kises  type  transformation  of  the 
governing  equations  to  a  stream  function  coordinate  system  and  an  expli¬ 
cit  finite-difference  solution  of  the  resulting  equations.  To  illustrate 
the  technique,  consider  the  axial  momentum  equation;  which,  for  axisym- 
metric  flow,  in  boundary  layer  form,  can  bo  written 
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Using  the  transformation 


if  dtp  -  p  u  r  dr 


(1-2) 


at  cons  tail  L  x,  Eq.  1-1  becomes 

3u  =  1  _3_  ( r,?_u_vTr.2  3u  I  _  ±  3j)  , 

3x  1 {)  ~di>  |  ij)  3Q  |  pu  3x  ' 

This  equation  is  valid  throughout  the  flnwfield  except  on  a  centerline 
where  both  r  and  if'  become  zero.  However,  r.he  limiting  form  of  Eq.  1-3 
for  the  centerline  can  be  obtained  through  use  of  L'Hospital's  rule, 
yielding 
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Similar  transformation  and  limiting  form  behavior  applies  to  the  other 
governing  equations  in  the  flow. 
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The  governing  partial  differential  equations  describing  the 
conservation  of  mass,  momentum,  and  energy  within  the  flowfield  and  along 
the  axis  can  be  readily  put  in  an  explicit  finite-difference  form.  A 
backward-difference  scheme  is  employed  for  trie  axial  derivatives  and  a 
central-difference  scheme  for  the  radial  derivatives. 


Consider  the  flowfield  divided 
into  a  gt'id  in  i|>,  x  coordinates. 
Then  the  derivatives  of  an  inde¬ 
pendent  variable,  say  F,  will  be 
evaluated  by 
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where 
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